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We report on our ongoing effort to understand quenched lattice QCD at finite baryon number density. The 
quenched theory is sensitive to the baryon mass both at strong coupling and in the scaling region. However, we 
find that the quenched model is pathological for fi > at fi = 6.0, in agreement with past Lanczos analyses 

of the Dirac operator. 



This note contains an upgrade of our work at 
f3 = 6.0 [1] to larger lattices, and a report of new 
results in the strong coupling limit. The consi- 
derations which motivated us to start a new pro- 
gram at finite density have been presented in de- 
tail in the past [1]. Our primary goal was to 
obtain accurate measurements of thermodynamic 
observables and spectra, which give complemen- 
tary information on the physics of the transition. 
We have had partial success, since the results 
were satisfactory up to /j, ~ m w /2. Particularly 
interesting were the results for the baryon mass. 
Its behaviour was found to be in agreement with 
simple potential models challenged by previous 
lattice studies. A very modest bending in the chi- 
ral condensate was observed for /j, ~ m w /2. Un- 
fortunately, in the interesting region, /j, > m w /2, 
the fluctuations inherent in the simulation were 
so large that it was impossible to obtain sensible 
results. Then, as a further step, we decided to 
investigate these problems by moving to bigger 
lattices, and by using a different source for the 
inversion of the Dirac operator. We have thus 
obtained results at [3 = 6.0 that we shall present 
below. 

We have also proposed an explanation for the 
problems afflicting the quenched model which re- 
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lies on the fact that poles in the fermion propa- 
gator for n ~ m w /2 can indeed be there on isola- 
ted configurations, as observed in past work [1], 
[2] , thus giving hints of deconfinement and chiral 
symmetry restoration. However, the contributi- 
ons of these poles should cancel in the ensemble 
average of physical quantities, if confinement is 
realized. In our new work, we have searched for 
direct evidence of such poles by building an un- 
physical operator, a 'baryonic pion', and we have 
found the pathologies expected for /j, > m w /2. 
These measurements were made at strong coup- 
ling, where, among other technical advantages, 
m w /2 and m^/3 are well separated, and the lat- 
tice data are also amenable to a cross check with 
analytic computations [3]. We begin now the ex- 
position of our results. 

At [3 = we have generated 20 configurations 
on a 8 3 x 16 lattice. We have inverted the Di- 
rac operator three times, for the three different 
boundary conditions [1], for each m, pt. The bare 
quark mass was .1, /j, ranged from 0. to 1.2, and 
the mean field prediction for the pseudocritical 
H is n ~ 0.6. We have measured the chiral con- 
densate, the number density, the energy density 
and the pion mass. The other masses are too 
heavy, and difficult to extract. Along with stan- 
dard observables, the unphysical operator G^Gj, 
(the 'funny' or 'baryonic' pion, as opposed to the 
real pion G^Gl^) was measured to obtain direct 
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Figure 1. Chiral condensate (circles, left scale) 
and number density (crosses, right scale) as a 
function of /j, at [3 = 0. The solid lines are to 
guide the eye. 
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Figure 2. 'Funny pion' propagator at [3 = 0. 

= 0., 1, .2, .3, .4, from top to bottom. Note the 
'baryon-like' behaviour at small /i, and the flat- 
ness for jj, = .4 



evidence for poles of the Dirac operator. 

We have observed a dramatic increase in the 
cpu time required for the inversion of the Di- 
rac operator when entering the /j, > m w /2 re- 
gion. This is consistent with the electrostatic ana- 
logy [2], since in this region the Dirac matrix is 
almost singular. For /j, > itjb/3 the number of 
iterations needed for the inversion is again very 
small. In this region we are inside the eigenvalue 
'crescent', and the first Brillouin zone is satura- 
ted: all the observables have their limiting va- 
lues. The chiral condensate is zero and the num- 
ber density is three (Fig.l). This tells us that the 
theory in the strong coupling limit 'knows' about 
the baryon mass. 

We note that the average value of the chiral 
condensate is consistent with the one previously 
reported in the literature. However, for /j, > m w /2 
the results for < ipip > obtained configuration by 
configuration spread over a wide range, which in- 
cludes both zero and the /j, = value. This sug- 
gests that we might have entered a metastability 
region rather than the symmetric phase, and that 
the identification of the onset for < ipip > with 
the pseudocritical point computed in a mean field 
approximation may not be correct. 

The behaviour of the 'funny' pion is very inte- 
resting. Its propagator has a baryon-like (i.e. T- 
asymmetric) behaviour for /j, < m w /2, and then 
flattens 'inside' the eigenvalue cloud (jjl > m w /2) 



(Fig. 2). We have clear evidence of zero modes 
in the quark propagator for /j, ~ m w /2, again in 
agreement with the eigenvalue picture. The ma- 
thematical and physical significance of this result 
is open to many interpretations which have been 
discussed and reviewed in the literature [1]. 

At [3 = 6.0 we have accumulated 30 indepen- 
dent configurations, on a 16 3 x 64 lattice. A sub- 
set of configurations was used for the inversion of 
the Dirac operator for each /j, value. The quark 
mass was .02. The role of the winding loops was 
monitored by doubling ( with respect to our pre- 
vious work) the temporal extent of the lattice. 
Also, we decided to use a wall-noisy source to 
reduce the non-positivity effects associated with 
the rigid wall we were using before. 

We confirm the decline observed in the chi- 
ral condensate (Fig. 3). The huge and unphysi- 
cal fluctuations in the pion propagator described 
in [1] disappear. However, the number of itera- 
tions required for the inversion when /j, exceeds 
mjr/2 remains big, of order 10 4 . This means that 
the Dirac operator is still nearly singular. The 
source affects the amplitudes, but not the poles, 
as expected. Despite the better behaviour of the 
propagators, the mass estimates were not com- 
pelling: it is well known that a point source, even 
with the noisy improvement, has poorer behavi- 
our than a rigid wall one. In particular, we were 
not able to measure the baryon mass. 



3 



T = 64 o. T = 33 □ 

CI . 1 rj 

; ■ — — T 

5£ o.io 

o. or, — 



Figure 3. Chiral condensate as a function of fi on 
the 16 3 x (32, 64) lattices, from a noisy estimator 
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Figure 4. Pion and nucleon masses as a function 
of fi on the 16 3 x (32, 64) lattices, from a (rigid, 
noisy) wall source. 



rations required for the inversion was still huge, 
meaning that we were still inside the eigenvalue 
cloud. In the free theory saturation occurs only 
for large /i, and we are probably observing a si- 
milar trend at [3 = 6.0, consistent with the obser- 
vations of [4]. So, the physics of the quenched 
model close to the continuum for /j, > m w /2 still 
evades us. 

In conclusion, is the finite density puzzle sol- 
ved? Is there a phase transition at /j, = mg/3, 
obscured by precursor phenomena? Or is the 
transition at /j, = m^/3 really preempted by a 
spurious one at m w /2? Clearly, something patho- 
logical is occuring at /j, ~ m^/2-it has prevented 
us from getting results for /j, > m w /2 at [3 = 6.0. 
But it is also clear that the theory 'knows' about 
the baryon mass, and that it behaves in a sensi- 
ble way at small values of the chemical potential. 
Our next step would be a full QCD simulation at 
finite temperature, and moderate values of chemi- 
cal potential, to be run in parallel with the quen- 
ched model. We look forward to phenomenologi- 
cally relevant results, and, by comparing the two 
models, to a better understanding of the failures 
of the quenched approximation. Of course, it is 
unclear whether simulation methods, designed for 
unquenched QCD, such as the complex Langevin 
algorithm, will work at nonzero /j, since they de- 
pend on the existence of the quark propagator 
configuration by configuration. 



We have some indication of an increase of the 
pion mass across the transition (Fig. 4), consistent 
with the decline of the chiral condensate, and sug- 
gestive of chiral symmetry restoration. The ef- 
fect is not dramatic ( The pion mass is basically 
constant within its statistical errors. Moreover, 
wall and noisy-wall sources give lower and up- 
per bounds to the pion mass, respectively, so one 
may be tempted to 'average' the results, thus ob- 
taining a constant pion mass), but it is corrobo- 
rated by the trend in the propagators. In fact, 
we have observed that the pion and the sigma 
propagators become nearly degenerate. 

Also in this case, as we did at [3 = 0., we tried 
to get inside the eigenvalue distribution by incre- 
asing jj, to jj, = 1.0. However, the number of ite- 
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